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The crossing from a static quark- ant iquark system to a 
system of two static-light mesons when the separation of the 
static quarks is increased is calculated in zero-temperature 
lattice QCD. The mixing of these two states is extracted from 
the lattice operators. We also discuss the breaking of an ex- 
cited string of a hybrid meson. 

PACS numbers: 12.38.Aw 12.38.Gc 13.25.-k 13.75.Lb 



I. INTRODUCTION 

The breaking of a long flux tube between two static 
quarks into a quark- ant iquark pair is one of the most 
fundamental phenomena in QCD. Because of its highly 
non-perturbative nature it has defied analytical calcula- 
tion, while its large scale, e.g. when compared to the sizes 
of composite particles in the theory, has caused difficul- 
ties in standard nonperturbative methods. Thus string 
breaking has remained a widely publicized feature of the 
strong interaction that has never, apart from rough mod- 
els, been reproduced from the theory. 

String breaking can occur in hadronic decays of QQ 
mesons and is especially relevant when this meson is lying 
close to a meson-antimeson (QqQq) threshold. For the 
heaviest quarks involved in these decays applying heavy 
quark effective theory is a reasonable approximation. 

Due to recent advances in both computational hard- 
ware and algorithms, much interest in lattice QCD has 
been devoted to attempts to observe string breaking. The 
direct approach of trying to see the flattening in the static 
QQ potential at large separation has been successful only 
at temperatures close to the critical one [1,2]. The failure 
of this Wilson loop method at zero temperature seems to 
be mainly due to the poor overlap of the operator (s) with 
the QqQq state [3,5]. In much more easily calculable ad- 
joint SU{2) and 5'[/(2)+Higgs models without fermions 
a variational approach with explicit inclusion of both the 
Wilson loop and scalar bound state operators has worked 
well [3-5]. In three-dimensional SU(2) with staggered 
fermions an improved action approach has been claimed 
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to be successful with just Wilson loops [6]. In QCD with 
fermions effective operators for QqQq systems are, how- 
ever, hard to implement; part of the problem is the ex- 
hausting computational effort required to get sufficient 
statistics for light quark propagators with conventional 
techniques for fermion matrix inversion. 

A new technique of calculating estimates of light quark 
propagators using Monte Carlo techniques on pseudo- 
fermionic field configurations [7] with maximal variance 
reduction [8] has been found to be very useful for systems 
including heavy quarks taken as static, such as single 
heavy-light mesons and baryons and also two heavy-light 
mesons [9,10]. The application of this method to the 
string breaking problem seems natural. 

Our previous work [9,10] has concentrated on bound 
states of two heavy-light mesons and mechanisms of their 
attraction for various values of light quark isospin Iq and 
spin Sq. Here we continue by studying the QqQq sys- 
tem for Iq^Sq = (0,1) together with the QQ system at 
distances around the string breaking point ri, ~ 1.2 fm, 
where the energies of the two systems are equal if there 
is no mixing. 



II. QUANTUM NUMBERS AND HYBRID 
MESONS 



When a quark and an antiquark are created from the 
vacuum they should have the 0++ quantum numbers of 
the vacuum. A quark-antiquark pair with J^^ = 0++ 
has lowest orbital angular momentum L = 1 and is in 
a spin triplet. This so-called Quark Pair Creation or 
^Pq model was combined with a harmonic oscillator flux 
tube model by Isgur et al. [11] to describe local breaking 
(formation) of a flux tube. In our calculation the sym- 
metries of the static approximation for the heavy quark 
automatically lead to only the light quark spin triplet 
being nonzero, which can be seen from the Dirac spin 
structure of the heavy-light diagram in Fig. 1. 

String breaking can also occur for hybrid mesons where 
the gluon field between two quarks is in an excited state. 
Table I presents the couplings of some low-lying glu- 
onic excitations to the quantum numbers of the resulting 
meson-antimeson in the static limit for the heavy quarks. 
In this limit CP and J^, where z is the interquark axis. 
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are conserved. The table lists the representations of these 
symmetries (with S,n, A corresponding to Jz = 0, 1, 2 
respectively and u corresponding to CP==bl). For 
the lowest-lying excitation, which has 11^^ symmetry with 
= 1, only non-zero angular momenta L -\- L' for the 
resulting mesons Bl', Bl would be allowed, as Sq has to 
be zero to generate a negative CP. 

Here, as for the other symmetries, /g = as we need 
the same flavour for the light quarks. The Iq = 1 cases 
do not correspond to spontaneous breaking of a string 
in vacuum but a correlation of a QQ-\-qq system with 
a Qq-\-Qq system at different time, i.e. the breaking of 
a string in the presence of a meson. When the table is 
extended to nonzero L for the heavy quark plus antiquark 
(as opposed to L > for a single meson) the P, C values 
get multiplied with (—1)^. The energy levels also change; 
these retardation effects have been found to be relatively 
smaU [12]. 

Previously string breaking in hybrid mesons has been 
discussed from a phenomenological point of view using an 
extension of the approach of Isgur a/., i.e. a nonrela- 
tivistic flux-tube model with decay operators from strong 
coupling limit of lattice gauge theory and heavy quark ex- 
pansion of QCD in Coulomb gauge [13]. From this model 
two selection rules were given, the first one agreeing with 
the Hu case in Table I; low- lying hybrids do not decay 
into identical mesons, the predominant channel being one 
s and one p-wave meson. The second rule prohibits de- 
cay of spin singlet states into only spin singlets, which is 
not relevant to our calculation as our heavy quark spin 
decouples. 

An important question for hybrid meson phenomenol- 
ogy is the nature of the lowest state for a given set of 
quantum numbers at a particular heavy quark separa- 
tion; a hybrid QQ meson, a ground-state QQ meson with 
a qq meson or a system of two heavy-light mesons. It is 
also useful to know the strength of mixing between these 
states. This information can be obtained, in principle, 
from lattice calculations and used to decide what sort 
of bound states are most likely to exist and what their 
decays will be (see also Ref. [14]). 



III. LATTICE CALCULATION 

We use SU(3) lattice QCD on a 16^ x 24 lattice with 
the Wilson gauge action and the Sheikholeslami-Wohlert 
quark action with a nonperturbative "clover coefficient" 
csw = 1-76 and (3 = b.2 with two degenerate flavours 
of both valence and sea quarks. The measurements were 
performed on 20 gauge configurations. The gauge con- 
figurations are the same as in Ref. [15] for = 0.1395. 
With these parameters we get a lattice spacing a ~ 0.14 
fm and meson mass ratio Mps/My = 0.72. 

Estimators of propagators of quarks from point n to 
point m can be obtained from pseudofermion fields 0. For 



each gauge configuration a sample of the pseudofermion 
fields is generated, and the propagators are then obtained 
by Monte Carlo integration [7] . Thus there is one Monte 
Carlo averaging for the gauge samples, and another one 
for the pseudofermion samples for each gauge sample. 
In order to reduce statistical variance of propagators a 
variance reduction method similar to multi-hit can be 
used [8]; such a reduction is essential in practice. Our 
variance reduction involves division of the lattice in two 
regions, whose boundary is kept fixed while the ^-fields 
inside are replaced by their multi-hit averages. We use 
24 pseudofermionic configurations per each gauge config- 
uration. 




FIG. 1. Diagrams involved in the calculation; the Wilson 
loop VK, the heavy- light correlator t/, the unconnected me- 
son- ant imeson D and the box diagram B. 

Figure 1 shows the diagrams involved in the calculation 
with the time axis in the horizontal direction. The solid 
lines are heavy quark propagators, which in the static 
approximation are just products of gauge field variables. 
The wiggly lines are light quark propagators, obtained 
essentially as a product of pseudofermionic variables from 
each end which have to be in different variance reduced 
regions [8]. 

In the large T limit both the quark- antiquark and two- 
meson operators should in principle approach 
with Eo{R) being the ground state of the system. In 
practice the Wilson loop has a very small overlap with 
the two-meson state, which leads to great practical diffi- 
culties in observing the flattening Eo{R) 2MQq from 
it at large R. The heavy-light term U is necessary to 
obtain the correct ground state by explicitly including 
both quark- antiquark and two-meson states, and allows 
us to measure their overlap, which is crucial for string 
breaking to happen. 

To estimate the ground (and excited) state energy of 
our observables we always use a variational basis formed 
from different degrees of spatial fuzzing of the operators. 
This allows the use of moderate values of T instead of 
the infinite time limit to reduce excited state contribu- 
tions. The resulting correlation matrix C{R^T) is then 
diagonalised to get the eigenenergies. 

Due to the variance reduction method dividing the lat- 
tice in two halves in the time direction, the box diagram 
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and the heavy-light correlator in Fig. 1 have to be turned 
"sideways" on the lattice; i.e., the time axis in the dia- 
grams is taken to be one of the spatial axes to keep the 
light quark propagators going from one variance reduced 
volume to another. This induces technical complications 
that greatly increase the memory and CPU demands of 
the measurement program. 

For two flavours the Iq = wavefunction is of the form 
(i^iZ+c^J) /a/2, which gives factors of 1, a/2, 2 for I), [/, B 
respectively. For light quark spin we get the triplet states 
as in Ref. [9]. 

In this first study we concentrate on the ground state 
breaking, i.e. the first row of Table I. Investigation of the 
hybrid meson breaking requires diagrams not included in 
Fig. 1, which involve the hybrid QQ and QQ + qq oper- 
ators. We estimate that for the 11^ excited state the ex- 
cited string breaking happens in the same distance range 
as for the ground state due to the non-zero momenta of 
the resulting mesons (masses taken from Ref. [8]), which 
makes it harder to obtain sufficient accuracy as the spa- 
tial operators for excitations involve subtractions rather 
than sums of lattice paths. 



IV. RESULTS 
A. Variational approach 

A full variational matrix involving the Wilson loop, 
heavy-light correlator and the QqQq correlators gives the 
ground and excited state energies and corresponding op- 
erator overlaps as a function of heavy quark separation, 
in analogue to the approach of Refs. [4,5] for the adjoint 
string breaking and Refs. [3] for the SU(2)+Higgs model. 
We use a local light quark creation (annihilation) opera- 
tor and an extended version where a fuzzed path of link 
variables with length two separates the operator from the 
heavy quark line. For the link variables involved in the 
QQ operators we have two fuzzing levels. The two QQ 
and three QqQq basis states then give a 5 x 5 correlation 
matrix C{R^T) that can be diagonalised. However, for 
our present statistics the full matrix gives a reasonable 
signal only for r < r^. 

In Figure 2 the results from a calculation using just the 
most fuzzed basis states for both QQ and QqQq (a 2 x 2 
matrix) are shown. At we would expect the ground 
and excited state energies to be separated by twice the 
mixing coefficient x (see below). We observe a larger 
separation which is presumably due to our statistics not 
being sufficient to give accurate plateaus for the energies. 

Although this full variational approach is in principle 
the most direct way to study string breaking, we find 
that it is possible to focus on string breaking explicitly, 
as we now discuss. 
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FIG. 2. Ground and excited state from a variational calcu- 
lation including QQ and QqQq operators. The highest fuzzed 
basis state for both is used here. The ground state of the 
Wilson loop and 2mQq are also shown. 



B. Mixing matrix element 

In full QCD, there is mixing of energy levels between 
states coupling to Wilson lines (flux tube) and QqQq 
states. To get the mixing matrix element the correla- 
tion between a Wilson line and a QqQq operator has 
to be considered. In order to study the operator mix- 
ing from this heavy-light correlator one needs to use re- 
sults (energies and couplings) from both diagonal opera- 
tors separately: thus from the Wilson loop (with ground 
state contribution given by W{T) = w"^ exp[—V{R)T]) 
and the unconnected QqQq correlator (eg. D{T) = 

exj[)[— M {R)T] from the ground state) where we use 
a variational basis to suppress excited states. 

The ground state contribution to the heavy-light cor- 
relator can then be written as 

T 

[/(T) =x(i?)^^e-^^^)^e-^^^)^^-^)c/+0(x3) (1) 
t=o 

In the quenched case the contributions from fermion 
loops inside the correlator are absent, removing the 0{x^) 
terms in Eq. 1. The box term is expressed in the same 
manner as 

ti=0t2>ti 

X e-^(^)(^-*^)(i + Oix^) (2) 

The operator mixing coefficient x for the QQ and 
QqQq states can be extracted from these expressions. 
Near the string breaking point (where V{R) = M{R)), 
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in the infinite time limit, only the ground state contribu- 
tions survive. We use 



U{T) 



fT/2 



^W{T)D{T) 1 + . . . + 



+ 0(x3) (3) 
+ 0(x^) (4) 



V D{T) ^i + ... + (T+l)/T 

The factors of / = exp(y(i?) - M{R)) account for de- 
partures from the string breaking point. 

In the quenched case there is no mixing between the 
energy levels of the quark- ant iquark and two-meson sys- 
tems, and X can be extracted using Eqs. 3,4. As x « 1 
the non-leading terms in the expressions for x are small 
and we may use also these formulas with our unquenched 
data - with a resulting decrease in errors compared to the 
full variational study of the preceeding subsection. 

Our assumption about neglecting excited state contri- 
butions can be tested by obtaining consistent results for x 
from both relations for several T values. To improve fur- 
ther on our estimate of x we diagonalise separately D 
and B to enhance the ground state contributions and use 
the first two diagonalisations to extract the ground state 
of U. Our results with bootstrap errors can be seen in 
table II. Assuming constant x for 0.99 fm < r < 1.31 
fm gives us a best estimate of x/a = 0.033(6)/a = 46(8) 
MeV. This is about half of the value of x = 100 MeV 
obtained using a strong coupling mixing model and the 
experimental T(4S') decay rate [16]. 

Our analysis shows that the string breaking matrix ele- 
ment is small but non-zero. We are able, however, to find 
two independent ways to estimate it (using all four dia- 
grams in Fig. 1) and we obtain x = 46(8) MeV with light 
quarks that are around the strange quark mass. This 
is the first non-perturbative determination from QCD of 
the string breaking matrix element. Because of its small 
value, direct observation of string breaking from the spec- 
trum is difficult to achieve. 
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TABLE I. Relation of gluonic excitations of a hybrid meson 
to ground-state properties of the meson pair resulting from 
string breaking, in the static limit for the heavy quarks. The 
last three columns refer to the quantum numbers of the light 
quarks in the meson-antimeson system. 



r (fm) 


Eq. 3 without O(x^) 


Eq. 4 without 0{x^) 


0.85 


0.040(3) 


0.081(8) 


0.99 


0.025(4) 


0.045(5) 


1.14 


0.025(3) 




1.31 


0.031(15) 


0.040(25) 



TABLE IL Operator mixing x extracted using a variational 
approach with two different formulas. The values are taken 
at T = 4 and have bootstrap errors. Physical dimensions are 
obtained by multiplying with = 1.39 GeV 
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